Accounting for Slow J /if} from B Decay 
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A slow J/ip excess exists in the inclusive B — > J/tp + X spectrum, and is indicative of some 
hadronic effect. From color octet nature of cc pair in b — > ccs decay, one such possibility would be 
B — > J/ip + K g decay, where K g is a hybrid resonance with sgq constituents. We show that a K g 
resonance of ~ 2 GeV mass and suitably broad width could be behind the excess. 

PACS numbers: 13.25.Hw, 14.40.Nd 
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Using 1.12 fb _1 data, the CLEO experiment published 
the inclusive B — ► J/ip + X spectrum 0. After sub- 
tracting XcU Xc2 and ip(2S) feeddown, there was a hint 
of excess events around Pj/^ ~ 0.5 GeV, where p*j/^ is 
J/ip momentum in T(4S f ) frame. With 6.2 fb _1 data, the 
Belle experiment presented |2| the inclusive B — ► J/ip+X 
spectrum. After feeddown subtraction, one could also 
infer Q an excess for V*j/^, ^ 0.8 GeV. Recently, the 
BaBar experiment published Q similar results based on 
20.3 fb _1 data, showing clear excess beyond NRQCD ex- 
pectations 0, of order 1CT 3 in rate, for V*jt^ ^ 0.8 GeV. 

As the excess involves slow moving J ftp mesons, it 
must have hadronic, rather than perturbative, origins. 
Various proposals have been advanced. The suggestion 
of _B — > J / ipAp J(| has been studied recently by BaBar ; 
the event rate at 10 -5 order cannot explain the slow J/ip 
excess. Intrinsic charm content of the B meson could 
lead to B — > J/ipD(+7r) final states (sj, which can in 
principle explain the data, but experimental studies are 
not yet forthcoming. If B — > J/ipDir dominates, the slow 
pion does not pair with D to form a D* , and would pose a 
challenge. Another possibility Q would be B — > J/tpK g , 
where K g is a hybrid meson with sgq constituents. A 
recent estimate Q suggests that the rate could be in the 
ballpark. 

In this note we take a heuristic approach to explore 
the last possibility. We find the hybrid scenario is indeed 
viable. We suggest the signature of B — ► J/ip + K + nn, 
where n cannot be more than a few, should be experi- 
mentally searched for. If the K+mr system tends to peak 
at some mass, but does not descend from D meson decay, 
then the hybrid meson picture could be substantiated. 

Let us visualize why a hybrid meson recoiling against a 
J/ip could be the right picture. In b — ► ccs decay, the cc 
pair is dominantly formed in a color octet configuration, 
hence charmonium production is color-suppressed. Imag- 
ine that, upon b quark weak decay, the c and c quarks 
are moving apart with more than ~ 1 GeV kinetic en- 
ergy. Soft "muck" effects cannot change the configura- 
tion, and the system would tend to break up into open 
charm meson pairs, resulting in D^D^ or DDK final 
states. But if the c and c momenta are relatively colin- 
ear, it can be viewed as a "proto-charmonium" . Because 



of the heaviness of m c , this small, dominantly color octet 
cc system would recoil against the s quark, again rela- 
tively unperturbed by the soft "muck" . By the time it 
separates from the s quark by order 1 fm, strong, non- 
perturbative effects set in: it has to hadronize. But since 
this is already a "proto-charmonium" , i.e. the spatial and 
spin wave-function already maps well onto some physi- 
cal charmonium state, the only problem is it has to shed 
color. An effective color octet charge is thus left to neu- 
tralize the sq system, and the simplest configuration is 
that of sgq, which we call an K g hybrid system. 

Whether a K g hybrid meson really exists becomes se- 
mantical. As visualized above, the leftover color octet 
charge with the color octet sq system is not by perturba- 
tive gluon emission @ , but by the fact that two separate 
color strings extend from the s and the q towards the 
point where the cc hadronization took place. It is plau- 
sible that this s-g-q string configuration could resonate 
and the amplitude gets enhanced, and energy-momentum 
is exchanged with the departing charmonium. If the K g 
resonance is of order 2 GeV 10 in mass, then the J/yb 
has to be slow. We have therefore constructed the phys- 
ical picture whereby slow J/ip (charmonium in general) 
can receive enhancement. 

Let us illustrate further this picture. The relevant ef- 
fective weak Hamiltonian is 

H w = -Lv c * s V cb (c 1 1 +c 2 2 ), 

O l{2) = 5„y(l- 75)c a(/3) 5^)7^(1 -75K- (1) 
The amplitude for B — > J/tpK decay is 

A(J/ipK) - ^V c * s V ch a 2 (J/iP\n»c\0){K\s lfi b\B), (2) 

where a 2 is the effective coefficient governing color- 
suppressed modes, which in the naive factorization limit 
is C2 + ci/3. Fitting data gives a 2 > 0.25, but factoriza- 
tion calculations tend to give lower values of a 2 < 0.2, 
suggesting that nonfactorized effects are important. 

We now infer from Eq. J3J) by analogy the formula for 
B — > J/ipK„. It is well known that 
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where T a is the color SU(3) generator. A color octet cc 
pair is favored. The "proto-J/V'" would still be produced 
by a cc vector current, but stripping off the extra color 
as it departs, one is left with a "constituent" gluon in 
association with the sq bilinear. The nonperturbative 
picture should be two (different colored) strings extend- 
ing from the point of departure of the J/tp towards the 
recoiling s and spectator q quarks. Using a factorization 
language for sake of illustration, the matrix element prod- 
uct, (J/ijg\(cT a c) v _ A \0) ((sq) s \{sT a b) v „ A \B), where g 
is a constituent gluon and (sq)g is a color octet quark 
pair, should be nonvanishing. The final state therefore 
has the J/tp recoiling against a color singlet sqg con- 
figuration, which could form a hybrid K g meson. We 
heuristically write the operator 2 (cT a c)v-A(s T a b)v-A 
as C7 M c [sgb]^, which we now take as our Ansatz. Thus, 
the formula for B — > J/ipK g becomes, 

A{J/i>K g ) cc ^V c * s V cbCl (J/^\c^c\0)(K g \[sgb}^\B), 

(4) 

where the [sgb] operator can convert the B meson into 
Kg, including matching its J p , in analogy with the sj^b 
current converting B into K(*>. Note that the decay am- 
plitude is proportional to Ci *~ 1 rather than a2, but 
there is some proportionality constant, hopefully of order 
1, from the Ansatz we made above. Due to the usual 
difficulty of hadronic physics and the model dependence 
that would necessarily arise, we do not attempt at cal- 
culating theoretically this proportionality constant, but 
turn to data for its determination. 

The simplest situation would be to have a K g meson 
with J p = T (note that there are no "exotic" kaon hy- 
brids). The contraction of vector current indices should 
then be very similar between Eqs. J2J) and 10} , and the 
B — ► J/ipK g (0 T ) decay rate is estimated to be 
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where kq is related to the aforementioned proportionality 
factor, but it also contains possible differences in B — > K g 
and B — > K form factors, and PK {g) is the momentum of 
K( g ) in the B decay frame. Since the hybrid K g meson 
is expected to be broad, the decay rate is modulated by 
the Breit-Wigner factor pl| 
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where q 2 = m 



i 2 B — IvtlbEji^ and Ej/^ is the J/tp 
energy in the B rest frame. Note that, to account for 
kinematic dependence a 



factor is used instead 



of m. Furthermore, 
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We now make a fit to the inclusive J/tp spectrum and 
see whether a single hybrid K g suffices to account for 
the observed excess. The direct (feeddown subtracted) 
B — > J/ipX data is taken from Ref. £|. BaBar uses 
NRQCD plus B -> J/tpK^ simulation results to fit their 
data. To simplify, we follow Refs. 0, @ , and use 



f(p)=N(p- Pmin)(p ~ Pmax) exp 



(p - pf 



, (8) 



to mimic the color-octet NRQCD and J/ipK^ 
components. We find (N, p min , p max , p, <r ) = 
(26, 0, 1.95, 1.21, 0.5) and (180, 1.2, 1.95, 1.65, 0.3), 
with all energy-momentum in GeV, give good accounts 
of the two contributions, respectively. 

Adding now the B — > J/tpK g contribution of Eq. 
we smear by a Gaussian with rms spread of 0.12 GeV 
to account for broadening in the T(4S') frame. Using 
|ci /a2 1 = 5 for illustration, and the isospin averaged 
B(B -> J/tpK) = 0.94 x 10" 3 , and allowing the NRQCD 
contribution to float in the new fit, we obtain 
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2.3, 



2.08 GeV, T ~ 72 MeV. 



(9) 



The NRQCD contribution is reduced by 12% with respect 
to Ref. U, the fitted B — * J/ipK g branching ratio is 
8.5 x 10~ 4 , and the spectrum is shown in Fig.^a). It is 
remarkable that a single O^ hybrid K g meson with mass 
- 2.1 GeV and width - 100 MeV could account for the 
observed slow J/tp excess. The fudge factor k ~ 2.3 
means our inference by analogy is in the right ballpark. 

Let us investigate the case for K g meson with J p = 1 T . 
Replacing K by K* in Eq. we estimate, in analogy 
to Eq. (0), the B — > J/iJjKg(l T ) decay rate to be 



4 



PKjL 
PK* 



|BW| 2 6(B-> J/^K*) : 



(10) 



K g and 



where K\ is analogous to kq, but now the B 
B — * K* form factor ratio can be rather complicated, be- 
cause of two possible helicity configurations. Absorbing 
all of this into «i, we retain linear momentum depen- 
dence corresponding to longitudinally polarized compo- 
nent, which is expected to be dominant from usual form 
factor models, as well as in the heavy quark limit. 

Performing a fit as before using isospin averaged 
B(B -> J/ipK*) = 1.35 x 10^3, we obtain 



m ~ 0.6, m ~ 2.05 GeV, T ~ 70 MeV. 



(11) 



where to, Tq are the mass and width of K g at q 2 



The NRQCD contribution is reduced by 11% with re- 
spect to Ref. 0, and the fitted B(B -> J/tpK g (l^)) is 
7.9 x 10 -4 . The spectrum is shown in Fig.^b), which is 
similar to Fig. QJa) since mass and width are almost the 
same. Note that the fudge factor k\ ~ 0.6 appears even 
more reasonable than the T case. 

Fig. ^ suggests that the fits may not yet be optimized 
for P*j/^ between 0.9-1.5 GeV. As we have allowed some 
freedom in the strength of the NRQCD contribution |{| , 
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FIG. 1: B — > J/ip + X decay spectrum. Data is from Ref. Q. 
The dotted and short-dashed lines correspond to J/ipK^ 
simulation Q and NRQCD component The long-dashed 
line corresponds to the B — > J/ipK g contribution for (a) 
and (b) case, with fitted mass and width given in text. 



we now allow its peak position to float as well. Fitting 
again, we find for case 



Kn ~ 2.0, m ~ 2.08 GeV, T ~ 147 MeV, 



(12) 



with the NRQCD contribution reduced by 13% with re- 
spect to Ref. [4( , and the parameter p shifted by 100 MeV 
to 1.31 GeV. The fitted B — > J/t/jK g branching ratio be- 
comes 12.9 x 10 -4 . For l^ case, we find 



Ki ~ 0.5, m ~ 2.03 GeV, T ~ 103 MeV, 



(13) 



with NRQCD contribution reduced by 10% with respect 
to Ref. 3, p shifted from 1.21 GeV to 1.29 GeV, and 
fitted B{B J/^K g {lT)) ~ 10.5 x 10~ 4 . The fitted 
spectrum is shown in Fig. [21 The remaining slight "dis- 
crepancy" can be attributed to the difference between a 
perturbative vs. hadronic approach, e.g. summing over 
J/tpKi 0, J/ipK%, etc. modes. 

We do not commit ourselves to what should be the 
lightest Kg hybrid state, or how large is the splitting 
for further excitations. It is gratifying, however, that 
the fitted masses of order 2-2.1 GeV is close to expec- 
tations The width of 70-150 MeV may seem nar- 
row, but kaonic hybrids have not been widely discussed 
in the literature, and the relative narrowness would make 
experimental identification easier. As for d ecay modes, 
we remark that flux tube models suggest [13| hybrids 
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FIG. 2: Same as Fig. 1 but allowing the NRQCD component 
to float for an improved fit. 



decay into final states with one excited meson. Thus, 
one should consider reconstructing K g in A'q 12 7t(p), 

•K' /o,i,2! maybe also Kfy final states. It would be fas- 
cinating if a heavy kaon resonance is found to be domi- 
nating the slow J ftp excess from B decay. 

We note that the K g width is larger for the improved 
fit of Fig. 2. We have checked that T ~ 250 MeV is 
possible, if the "peak position" parameter p is allowed to 
shift slightly higher, to 1.34 GeV and 1.31 GeV, respec- 
tively, for the T and 1 T cases. Thus, the narrowness of 
r ~ 70 MeV of Eqs. 10 and (JTTJ may be an artefact of 
trying to mimic the NRQCD result of Ref. Q . The latter 
work subtracted B — > J/ipK and J/ipK*, an approach 
BaBar adopted, but it was done before the Belle mea- 
surement [l3 of B(B -> (1270)) ~ 1.55 x 10~ 3 
(isospin averaged), which is comparable to B — > J/ipK 
and J/tJjK*. An update of Ref. 5| would be helpful. 

The 0~ and 1 — hybrid quantum numbers allow us to 
make some insight into the possible cause of sizable non- 
factorizable contributions to B — > J/ipK^*\ J/tpKi de- 
cay. The physical state may have a hybrid Fock 
component, \K) — c$\K^) + sg\Kg). One then finds 
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where a 2 ac ' is from factorization calculations, typically of 
order 0.15-0.2. One sees from our fitted n values that 
a hybrid admixture of a few % to no more than 10% 
in the K 1 -*' wavefunction can suffice to account for the 



4 



large ~ 0.25-0.3 extracted from data, in large part 
because of gaining the Cx/a^ factor. 

We comment on a recent proposal that slow J Jib from 
B decay could arise from cqqc four quark states [1J| . The 
physical picture would be that the color octet cc picks up 
a color octet qq pair as it hadronizes. Ref. [3 gave some 
arguments for why J jib may end up slow. It may not be 
so easy to distinguish this proposal from the present one, 
as both lead to J jib + K + nir final states, and in fact 
both mechanisms may well be at work concurrently. To 
distinguish the two mechanisms, one would have to check 
whether the charmonium system is resonating with some 
of the recoil hadrons. For K g mechanism, it would be 
helpful if Kq 1 2 ^(p)j K^fo,i,2, Kf% dominance in final 
state is borne out. We remark that for most pictures, 



one should also find excess in slow r\ c mesons. 

In summary, we have illustrated that a single hybrid 
Kg state recoiling against a J/tp could explain the slow 
J jib excess observed in B decay with rate at 10~ 3 order. 
The fitted K g mass is of order 2.05-2.1 GeV, with width 
of order 70-150 MeV, but could be broader. Experimen- 
tal signature would be to reconstruct K g — > K + rnr, 
probably in Kq 1 2 ^(p), K^* Jq,i,2, Kf? configurations. 
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